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In this paper, a new gain-scheduling control design 
approach is proposed by combining LPV control the- 
ory with interpolation techniques. The improvement 
of gain-scheduled controllers can be achieved from lo- 
cal synthesis of Lyapunov functions and continuous 
construction of a global Lyapunov function by interpo- 
lation. It has been shown that this combined LPV con- 
trol design scheme is capable of improving closed-loop 
performance derived from local performance improve- 
ment. The gain of the LPV controller will also change 
continuously across parameter space. The advantages 
of the newly proposed LPV control is demonstrated 
through a detailed AMB controller design example. 

Introduction 

The gain-scheduling approach is perhaps one of the 
most popular nonlinear control design techniques that 
has been widely used in fields ranging from aerospace 
to process control. Although it seems to work well in 
practice, this heuristic design procedure does not take 
the parameter variations into account. 1,2 In its early 
practice, the control design came with virtually no 
guarantee on performance, robustness, or even nom- 
inal stability. Recently, a systematic gain-scheduling 
design technique was developed in the form of linear 
parameter-varying (LPV) control theory. 3 " 8 This class 
of systems is different from the standard linear time- 
varying counterpart due to the causal dependence of 
its controller gains on the variations of the plant dy- 
namics. The implications of parameter-dependent sys- 
tems theory for gain scheduling is obvious because gain 
scheduling conceptually involves a linear, parameter- 
dependent plant. The LPV design technique provides 
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guaranteed stability and performance properties, and 
simplifies the interpolation and realization problems 
associated with conventional gain-scheduling meth- 
ods. However, one potential problem associated with 
parameter-dependent Lyapunov function approach is 
the lack of guidance for choosing “right” basis func- 
tions to parameterize infinite dimensional functional 
space. 

Interpolation is an important step toward synthesis 
of gain-scheduled controllers and has not been ade- 
quately addressed in a systematic way. Some ad-hoc 
interpolation techniques have been proposed in the 
past: (1) linear interpolation of poles, zeros, and gains 
of local controllers; (2) linear interpolation of solu- 
tions of Riccati equations; and (3) linear interpolation 
of state-space matrices of balanced controller real- 
izations. These approaches are intuitively appealing 
but could generate destabilizing controllers. Stilwell 
and Rugh 9 proposed a gain-scheduled control design 
based on interpolation techniques. The interpolated 
parameter-varying controller preserves point- wise sta- 
bility of local LTI controllers for all frozen parameter 
values. However, the global control does not provide a 
priori stability and performance guarantee from locally 
designed controllers for fast time-varying parameters. 

In this research, we propose an interpolating LPV 
control approach by combining LPV control theory 
with interpolation techniques. The interpolated LPV 
control design is capable of improving controlled per- 
formance by finding the most appropriate Lyapunov 
functions in a local sense. Moreover, the local stabil- 
ity property will be extended to the entire parameter 
range by a globally constructed Lyapunov function 
through interpolation. The proposed interpolating 
LPV control technique has great potential to many- 
industrial applications including active magnetic bear- 
ings (AMBs), for which the controller adjusts its gain 
based on changing rotor speed and provides accurate 
control of flexible mode variation in terms of natural 
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frequency, damping, and mode shape. In this research, 
we use a simple rigid magnetic bearing example to 
demonstrate advantages over existing gain-scheduling 
approaches. 

Improved LPV Analysis Condition 

Consider an LPV system 


X (t) 


\Aipit)) 

BiPim 

~xit) 

fit). 


[C(p(t)) 

Dipit)) \ 

d (t ). 


matrix function Qx{p) > 0, if there exists a group 
of continuously differentiable positive-definite matrix 
functions Xi(p),i = 1,2 such that for any p £ V % , 


-(A: i fp)xfip) + x l (p)Mp)\ 

{ ±v (Qx(j>) + J 

B T ( P )X i (p) 

C(p) 

and for all p € V 12 


Xiip)B(p) 

c T { P ) 

- 7 il 

D T (p) 

D{p) 

~lil 
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where x,x € R", d € R" d , e € R nc . All 
the matrices have compatible dimensions. It is as- 
sumed that the vector-valued parameter p evolves 
continuously over time and its range is limited to 
a compact subset V C R“. Its time derivative is 
bounded and satisfies the constraint —Vi < pi < 
vi ,i = 1,2, ,s. For notational purposes, V := 

{v : -Vi<Vi<Vi,i = l,2,---,s}, where V is a given 
convex polytope in R s that contains the origin. Given 
the sets V and V, we define the parameter p- variation 
set as 

Tf, = {pe C 1 (R+. R 5 ) : pit) 6 V, pit) G V, Vi > 0} 

( 2 ) 

Therefore the dynamics of the LPV system are char- 
acterized by the parameter value p and its variation 
along time. Previous research on LPV control the- 
ory mainly focused on a single Lyapunov function 
(quadratic or parameter-dependent) over the entire 
parameter set. For the given LPV system, it is clear 
that the achievable performance relies on the choice of 
the Lyapunov function. However, it would be bene- 
ficial if one can analyze the performance of the LPV 
system over different parameter ranges using different 
Lyapunov functions with a stability guarantee. 

To simplify the presentation, let us assume that the 
parameter set has dimension one. That is, V C R 1 . 
Suppose V \ , V 2 is an overlapped partition of the pa- 
rameter set V, and define the intersection of V 1 and 
V 2 as 

v 12 = v l nv 2 ( 3 ) 

Then the parameter space V is the union of two subsets 
V 1 ,V 2 . Through linear interpolation of the Lyapunov 
functions Xi ip), X\ ip) over subsets V 1 and V 2 , one 
would obtain a continuous Lyapunov function over 
the entire parameter space. The continuity property 
of the Lyapunov function is important to guarantee 
monotonic decrease of its value over any allowable tra- 
jectories dictated by the LPV system dynamics. A 
new LPV stability and performance analysis condition 
based on the partitioned parameter subsets is then pro- 
posed. 

Theorem 1 For the parameter set V = [p, p] with 
■its overlapped partition V 1 = [p, 6] and V 2 = [a, p] 
( a < b). the LPV system (1), and a given symmetric 


X 2 (p) - A f(p) + (b- a)Qx(p) > 0 (5) 

X 2 (p) ~ Xi ip) - ib - a)Qxip) < 0 (6) 

then the LPV system (1) is exponentially stable and 
its induced C 2 norm with x(0) = 0 is bounded by 7 = 
max{7i , 72 } • 


Given the above conditions over each parameter sub- 
set, a continuous Lyapunov function in the form of 


V(x,p) = x 1 


b p \x 1 ip) + (?— 7 )X 2 (p) 


b - a 


b - a 


(7) 

can be constructed and used to verify the desired sta- 
bility and performance properties of the system for any 
parameter trajectories within Tf,. Then a non-smooth 
dissipative system theory is applied to address possible 
discontinuity of the Lyapunov function derivative at 
the boundary of the parameter subsets. 10 The perfor- 
mance bound derived in Theorem 1 states the “worst- 
case” performance, which is not often achievable. For 
example, if the parameter trajectory is constrained in 
the subset V 1 , then the system’s performance is over- 
bounded by 7 i, which could be less than max {71, 72}. 

Theoretically, this theorem simply states the fact 
that the performance level of the LPV system relates 
to the existence of a continuous Lyapunov function. 
However, in practical situations, it is usually hard 
to identify a suitable parameter-dependent Lyapunov 
function without resorting to a global search. In this 
sense, the improved analysis condition will be help- 
ful to find a sharper performance bound through local 
studies. Then the global Lyapunov function will be 
constructed by interpolation. 


Interpolating LPV Control Synthesis 

Given a generalized open-loop LPV system 


x(t)' 


-4 (p(t)) Biipit)) B 2 (p(t)) 

x(t) 

e(t) 

= 

Ci(p(t)) Dnipit)) D 12 (p(t)) 

d(t ) 

U(oJ 


C 2 (p(t)) D 2 l {p(t)) D 22 {p{t))_ 
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( 8 ) 


where all of the state-space data are continuous func- 
tions of the scheduling parameter p, it is assumed that 
the parameter trajectory resides in the set Tf,. For 
simplicity, we assume that 
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(Al) (A(p), Bffp), C 2 (p)) is parameter-dependently 
stabilizable and detectable for all p, 

(A2) D 12 {p) and D 2 ffp) have full column and row 
rank respectively for all p , 

(A3) Dn(p) = 0 and D 22 (p) = 0. 


Theorem 2 For a partition of parameter space V = 
[p, p] as subsets V 1 = [p, 6] and V 2 = [a, p] (a < b), 
and given Qr,Qs > 0 matrix functions, if one of the 
following equivalent conditions are satisfied: 

1. there exist continuously differentiable matrix func- 
tions Rffp), Sffp), i = 1,2 such that for p £ V 1 


We consider the class of interpolating LPV con- 
troller in the form of 


X k (t) 


'A k (p{t),p(t)) 

B k (p(t),p(t)) 

Xk(t) 

_u(t) _ 


C k (p(t),p(t)) 

D k (p(t),p(t)) 

[ y(t)\ 


(9) 


In general, the controller gain is a continuous function 
of parameter p and its derivative. The control de- 
sign objective is to minimize the energy (£ 2 norm) of 
the output signal e(f) of the closed-loop system in the 
presence of bounded energy disturbance d(t). Specif- 
ically for the magnetic bearing control problem, we 
would like to synthesize a gain-scheduled control with 
combined disturbance rejection, gyroscopic compensa- 
tion, and automatic balancing capability. This can be 
formulated as an optimization problem with minimiza- 
tion of displacement of the rotor from its centerline 
at selected points subject to unknown torque distur- 
bances. 

Define x ^ := [ x T xj], then the closed-loop system 
can be written as 


Xd 


A c i (p, p) 
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Xd 
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Cd(pff) 
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c\ 
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Ak Bk 

0 / 

0 


C k D k 

c 2 o 
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( 11 ) 


Next, we propose a synthesis condition for an 
improved LPV controller using multiple parameter- 
dependent Lyapunov functions and an interpolation 
scheme. For clarity, only one parameter is considered 
and the parameter space V is covered by two over- 
lapped subsets V 1 and V 2 . For each parameter subset, 
we seek to design one LPV controller as stated in the 
form of eqn. (9). The overall gain-scheduling con- 
troller is then constructed by interpolating local LPV 
controllers. It is clear that the global controller is ca- 
pable of achieving tighter performance due to smaller 
parameter range. However, a critical issue associated 
with the proposed controller interpolation scheme is 
the stability of the global LPV controller. This will 
be guaranteed by constructing a globally continuous 
Lyapunov function over the entire parameter set. 



with N R {p) = Ker [Bj (p) 


DUp) 0 \,N s {p) = 


Ker [C 2 {p) D 2 i(p) 0], 

and 

for p £ V 12 
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(15) 

R2(P) 

-Rffp) - 

- (b - 
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s 2 (p) 

+ 

co 

1 

(b- 

a)Qs > 0 

(17) 

S 2 (p) 

-Sffp)- 

(b~ 
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(18) 


2. there exist continuously differentiable matrix func- 
tions Rffp), Sffp ) > 0 ,i = 1,2 and continuous 
matrix functions A^ff ,B k \C k \i = 1,2 such that 
for p £ V 1 , 



SjBi + B ( k i} D 21 
B l 

-iff 
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cl 

RiC r + C ( ff )T Df 2 
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-iff 


< 0 


'Ri(p) I 
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and for p £ V 12 


> 0 
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( 20 ) 


Rffp) - Rffp) + (b-a)Q R >Q (21) 
R 2 (p) - Rffp) - (b - a)Q R <0 (22) 

Sffp) - S, (p) + (b - o)Qs > 0 (23) 

Sffp)-Sffp)-{b-a)Qs< 0, (24) 
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then the closed-loop LPV system can be stabilized by 
a continuous LPV controller with induced C 2 perfor- 
mance less than 7 = max{7 1 ,7 2 }. 

Furthermore, let 


7 (P) 


7i 

(&)l. + (S )-« 


72 


p < a 
a < p <b 
p>b 


R(P),S( P ) 


(25) 

Ri(p),Si(p) p<a 

v / / \ a < p <1 

+ (^){R 2 (p),s- 2 (p)} 

R2(p),S 2 (p) p>b 

(26) 


and M(p)N T (p) := I — R(p)S(p), calculate 

A T (p) 


j , .. , dR . dM T 

Mp,p) = s(p)- + Mp)— 


(27) 


Cf(p)] 

x [ B i T W . 

4(P) = - [7(p)C 2 T (p) + S^B^p)^^)] 

x [2? 21 (p)DT(p)]- 1 (28) 

Cfc(p) =? - [-D^(p)-D 12 (p)] _1 

x [j(p)Bj (p) + ^(pJCj (p)fl(p)] (29) 

Dk(p)=0 (30) 


Then the interpolated LPV controller K-p will be con- 
structed as 


Ak(p,P ) = N-'ip) {i*(p,p) - S(p)B 2 (p)C k (p) 


-B k {p)C 2 {p)R{p) - 5(p)A(p)7?(p)} M- 

" T (P ) 

(31) 

B k (p)=N- 1 (p)B k (p ) 

(32) 

C k (p)=C k (p)Ar T (p) 

(33) 

D k (p)= 0. 

(34) 


It is noted that we ask for different performance 
levels over each parameter subregion. The proposed 
interpolation scheme could provide a global stabiliz- 
ing LPV controller with the potential to improve its 
local performance. The global Lyapunov function for 
the closed-loop LPV system is derived from the matrix 
functions Ri(p) and Sj(p). Because of the interpola- 
tion scheme used, it relaxes continuity requirement of 
individual Lyapunov functions over the intersections 
of parameter subsets. In particular, only finite vari- 
ation of Lyapunov functions over intersected regions 


is necessary and the amount of allowable change de- 
pends on the maximal parameter variation. When 
v — » 0 (slowly varying case), the conditions (16)-(18) 
will have no effect on the synthesis result. Thus the 
synthesis condition over each subset is essentially de- 
coupled and can be solved independently. 

The solvability conditions in eqns. (16)-(18) 
are clearly infinite-dimensional, as is the solution 
space. To approximate, we restrict the search of the 
parameter-dependent Lyapunov function to a span of 
finite numbers of basis functions. That is, let 

N f Ng 

Ri(p ) = E Si{p) = E 9f(p)Sf (35) 

k=l k= 1 

where { //"' (ft ) } ][' , and {sf (p)}E are user-specified 
scalar basis functions. Note that different basis func- 
tions may be used over different parameter subsets. 
Rlsf are new optimization variables to be deter- 
mined. After such a parameterization, the LPV 
synthesis conditions can be solved using a gridding 
method over each parameter subset. 

Advanced AMB Gain- Scheduling 
Control 

Active magnetic bearings (AMBs) use an electro- 
magnetic force to provide noncontact support for ro- 
tors in high-speed rotating machinery. AMBs have 
several unique characteristics that make them well 
suited for aerospace applications: they provide non- 
contacting support of the flywheel rotor, virtually 
eliminating bearing friction and wear, and they elim- 
inate concerns with lubrication, heating, and power 
consumption typical of standard bearing systems. Ad- 
ditionally, AMBs can effectively eliminate synchronous 
vibrations associated with mass imbalance and shaft 
run out, making them highly desirable for spacecraft 
applications where pointing accuracy is critical. Mag- 
netic bearings are well suited as flywheels in replace- 
ment of chemical batteries to store energy on a space- 
craft, 11 and they provide integrated attitude control 
and momentum management functionality. 12-14 Al- 
though this technique has tremendous potential for a 
variety of industrial applications, AMBs are open-loop 
unstable, thus making the controller design problem a 
challenge. Moreover, the flexibility of high-speed ro- 
tors adds to the complexity of control design task. 

Most controllers in use today for magnetic bear- 
ings were designed using PID strategies. However, 
it is difficult to satisfy the stringent performance re- 
quirements with PID control. Mohamed and Busch- 
Vishmiac 15 and Matsumura et al. 16 designed gain- 
scheduled TLoo controllers for rigid magnetic bear- 
ings utilizing the stabilizing controller parameteriza- 
tion with the free parameter Q playing the schedul- 
ing mechanisms. Unfortunately, this ad-hoc gain- 
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scheduling design technique suffers from slow vary- 
ing parameter requirements. The resulting gain- 
scheduling controller could render the closed-loop sys- 
tem unstable when the rotor speed changes rapidly. 16 
Rotor flexibility was considered using the //-synthesis 
technique. 17 It was shown that the p controller exhib- 
ited much greater stability robustness to variations in 
rotor mass. However, the controller caused significant 
performance degradation when the actual rotor speed 
differs from the pre-specified design condition. One 
way to address this problem is to design a series of con- 
trollers for each operating speed and then interpolate 
between these controllers. 16 However, this approach 
does not provide any stability and performance guar- 
antee. In Tsiotras and Mason 18 and Tsiotras and 
Knospe, 19 the LPV control theory was applied to ad- 
dress variable rotor operating speeds using a rigid rotor 
model and conservatism of LPV control design was 
recognized. Similar controllers were also developed by 
Sivrioglu and Nonami. 20 

Rotor Dynamics Modeling 

Owing to the linear dependence of the rotor speed 
in the plant dynamics, the nonlinear gyroscopic equa- 
tions of AMB can be simplified to a set of linear, 
time- varying differential equations. However, the ro- 
tor dynamics are inherently unstable; that is, even 
small unbalanced masses can create large synchronous 
disturbances with the same frequency as the rotor’s. 
Therefore it is necessary to develop a gain-scheduled 
controller capable of rejecting the periodic disturbance 
as the rotor speed changes. 

The rotational motion of a magnetic bearing can be 
derived from its rigid body dynamics, which is 15, 18 

9 = + ~r{fri - fr 2 + fn ~ fa + fdo) (36) 

t/f fJ y 

$ = t-y-9 + -r(fr3 - fr4 + fu ~ ft 3 + fdxb) (37) 


Parameter 

Value 

A area of each pole 

1531.79mm 2 

h pole width 

40.00mm 

Gq nominal gap 

0.55 mm 

J r radial moment of inertia 

0.333 kg ■ m 2 

J a axial moment of inertia 

0.0136 kg ■ m 2 

£ half the length of the shaft 

0.13 m 

k 

4.6755576 x 10® 

N number of coil turns 

400 

R coil resistance 

14.6 Ohm 

$o nominal airgap 

2.09 x 10 ~ 4 Wb 


Table 1 Magnetic bearing parameters. 


In the absence of disturbances and modeling errors, 
the above equation specifies an equilibrium. Lineariz- 
ing the nonlinear equations at the equilibrium, we 
obtain 

0 = + -f (~ 4c -^ + 2 Cl 0 e + f de ) (40) 

•J y i/f 

fp = ^+^(-4c 2 ^ + 2ci^ + /^) (41) 

N(j>e = eg + 2 d 2 £0 - di4>e (42) 

Nfa = + 2do(ip - di(j>^ (43) 


where <j>e,<f> $ are the differential magnetic flux from 
electromagnetic pairs, and e$,e$ are the correspond- 
ing differences of electric voltage. The constants 
Ci , Co , di , d ,2 depend on <S> 0 , Gq, R, A, N, v 0 and the ge- 
ometry of the bearings as follows 


ci = 2fc4>o 



2fc$g 
° 2 7 rh 


di 


2RG 0 2R% J r 

u 0 AN ’ d 2 ~ u 0 AN ’ m ~ £ ' 


(44) 

(45) 


The imbalance forces fee and fdxi> are typically mod- 
eled as sinusoidal disturbances, and are given by 


where 6, ip are the Euler angles denoting the orienta- 
tion of the rotor centerline. J a , J r are the moment of 
inertia of the rotor in axial and radial directions, re- 
spectively. The parameter p denotes the rotor speed. 
The magnetic forces generated by four pairs of elec- 
tromagnets are denoted by f r i,fn for i = 1,2, 3, 4. 
fdd i fdip are disturbance forces caused by gravity, mod- 
eling errors, imbalances, etc. 

The electromagnetic force fj is related to the voltage 
Cj across the jth coil through the magnetic flux tpj by 
the equations 


fj = H' 


. 2 g: 

1 + 

7TA 

2 R 


e,- = -| f—g j<p 

dt v 0 Al\ JJ ' 


(38) 

(39) 


The following system parameters are chosen for the 
active magnetic bearing example (Table 1). 


he = ( Jr ^ Ja \ J 2 Tcos(pt), (46) 

fdj, = ( ' Jr ( Ja p 2 rsm(pt). (47) 


In automatic balancing design, the imbalance forces 
will be treated as a sinusoidal sensor noise n T = 
[£t sin(pt + A) £rcos(pt + A)] on the measured rotor 
displacement, which can be written in a state-space 
form of 


%di — 


a 

P 

T 

0 


~P 

a 


%di + 



O' 

1 


%di 


(48) 

(49) 


where a = —0.001. 

Let x T = j (0 tip td tip 4>e (p^ xai , w = d 
and u T = [eg e^\. Combining rotor dynamics and the 
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disturbance equation, the linearized equation can be 
written as 



x = A(p) x + Biw + B 2 u 


(50) 


z = C\x + Dnw + -D 12 M 


(51) 


y = C 2 x + D 21 w + D 22 u 


(52) 

■ 0 

0 1 

0 

0 

0 

0 0 ' 

0 

0 0 

1 

0 

0 

0 0 

_ 4c2 

0 0 

_£Ja 

Jr 

0 

2ci 

0 

0 0 

o" 

_ 4C2 pjg 

d 

2ci 

0 0 

2d? 

N 

0 " d 

0 

_iti 

N 

d 

0 0 

0 

2d? (A 

N u 

0 

0 

di 

N 

0 0 

0 

0 0 

0 

0 

0 

a -p 

0 

0 0 

0 

0 

0 

P o - . 






(53) 

X 

to 

O 

’ S2 = ^ 

X 

0 




1 

0 

h 

_02x2_ 



(54) 

h O 2 X 6 n 

0 2 x8 ’ DU - 

X 

0 

D\2 

02x2 

h 



(55) 


C -2 — \l -2 02x4 h] ; D 2 I = O2X 1 , D 22 — 02x2 

(56) 

Note that the linearized rotor dynamic equation is in 
the form of a linear parameter-varying system with 
the rotor speed serving as the parameter. For gain- 
scheduled control, the rotational speed is assumed to 
be available in real-time for controller use. 



Fig. 1 Weighted open-loop interconnection for the 
magnetic bearing system. 

range. The LPV synthesis problem can be solved us- 
ing either a single or parameter-dependent quadratic 
Lyapunov function over all gridding points in a one- 
dimensional parameter space. 5,8 The performance ob- 
tained through a single quadratic Lyapunov function 
(SQLF) approach is 7 = 320.78, which is quite con- 
servative. Also, multiple LPV control syntheses are 
conducted with the whole parameter space divided 
into two subsets ([315, 720] U [700, 1100]) and four sub- 
sets ([315,520] U [500,720] U [700,920] U [900,1100]), 
respectively. Then five points are used to grid each 
parameter subspace uniformly. The free parameters 
Qn and Qs are chosen as O.I/12 for p e V 1 C 1 "P J or 
zero otherwise. In each partitioned parameter space 
case, two sets of identical basis functions are used to 
parameterize the functional space, 


Interpolating LPV Control Design 

The design objectives of the LPV controller are 
to asymptotically stabilize the system over the whole 
range of rotor speeds and to minimize an error signal 
representing a weighted sum of the forces at the bear- 
ings, the gap displacement at the bearing locations, 
and the control input used. 

Although the scheduling variable is time-varying in 
the LPV dynamics, it is simply treated as a fixed pa- 
rameter in the design stage. The design objective is 
quantified from a frozen parameter design viewpoint 
by weighting functions, and the weighted open-loop 
interconnection is given in Figure 1. 

In Figure 1 the weighting functions are chosen as 


W e (s) 

W u (s) 


200(s + 100) 
s + 0.01 2 ’ 

0.001s 
0.05s + 1 2 ’ 


W n (s) =0.001 


(57) 

(58) 

(59) 


The rotor speed is assumed to change freely be- 
tween 315rod/s to 1100 rad/s. The rotor dynam- 
ics exhibit significant gyroscopic effects in this speed 


flip) = 1, (60) 

9i(p) = h 9i(p)=P (61) 

for each i. Note that the basis function for Ri(p) 
is chosen as a constant over each parameter subset. 
This will only cause small performance degradation 
compared with more complicated basis function se- 
lections. Since both R and S become functions of 
scheduling parameters after interpolation, the inter- 
polated LPV controller gain will depend on both the 
parameter and its derivative. The performance level 
of the interpolating LPV controller versus number of 
parameter subsets under different parameter variation 
rates is shown in Table 2. Note that the 7 value 
represents the “worst-case” interpolating LPV control 
performance over different parameter subsets. The ac- 
tual performance over each parameter subset could be 
less than the “worst-case” performance. 

For comparison, a single parameter-dependent Lya- 
punov function is also considered to demonstrate the 
sharpness of the newly proposed LPV synthesis condi- 
tion. The calculated performance bound is about the 
same as the “worst-case” performance of interpolated 
LPV cases. 
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Variation rate 

PDLF 

Interpolating LPV 
2 subsets 4 subsets 

1 

9.863 

9.936 

9.954 

100 

10.064 

10.107 

10.149 

10000 

141.24 

144.43 

143.50 


Table 2 Achievable induced-/^ performance us- 
ing single and multiple parameter-dependent Lya- 
punov functions (PDLFs). 


Simulation Results 

Next, the LPV controller synthesized using two sub- 
sets is used for comparison study and simulation work. 

To compare the LPV controller with the optimal 
H oo controller, eight points are chosen from the param- 
eter set V and H oo controller is synthesized for each 
operating condition. Then their Hoo norm values are 
compared with the sub-optimal control performance 
that is achieved by LPV controllers evaluated at fixed 
parameter values. These are provided in the Table 3. 


p(rad/s) 

' Hoo 

Interpolating LPV 

315 

9.778 

9.833 

427.1 

7.487 

8.147 

539.3 

6.149 

6.719 

651.4 

5.289 

5.772 

763.6 

4.696 

4.988 

875.7 

4.264 

4.831 

987.9 

3.936 

4.600 

1100 

3.678 

4.383 






Table 3 Frozen optimal/LPV closed- loop "Hex, 
norm. 

Since the LPV controller is designed for a range of 
parameters, it is not surprising that they are only sub- 
optimal for each fixed parameter value. However, it 
is observed that the optimal performance level using 
Tfoo is very close to the achievable performance de- 
rived from switching the LPV controller at each frozen 
parameter. However, the optimal controller is highly 
tuned to its designated rotor speed. When a magnetic 
bearing operates at a different rotational speed, signif- 
icant performance degradation or even loss of stability 
is expected. For example, the optimal Hoo controller 
designed for p = 700 rad,/ s results in a performance 
level 7 = 1.16 x 10 6 when the rotor is actually rotat- 
ing with the speed of 1100 rad/s. 

Finally, we compare the closed-loop step response 
for both optimal and LPV controllers. The time re- 
sponses for the optimal controller and the LPV con- 
troller are shown in Figures 2 and Figure 3, respec- 
tively. It can be seen that the step responses of rotor 
displacement from both controllers are quite similar, 
whereas the control actions are slightly different due 
to sub-optimality of LPV controllers. 

After analyzing the frozen point LPV controller 
property, we are ready to do some nonlinear simula- 


(b) control force U \ , u -2 

Fig. 2 Fixed parameter unit-step response with 
Haa optimal controller. 

tions for AMB with time-varying rotational speed. A 
time-varying rotor speed profile is chosen as 

f sin(680f) 0 < t < 0.5s 

d{t) = < sin(u)t), to = 655 + 50f 0.5s < t < 1.5s 
[ sin(730 1) 1.5s < t < 2s 

(62) 

Note that the rotor speed trajectory is deliber- 
ately chosen to cross the intersection of two subsets 
[315,720] and [700,1100] of the scheduling parame- 
ter to illustrate LPV controller interpolating effect. 
Disturbance d = 1.3 x 10~° is used for the simula- 
tion purpose. Except for small glitches during the 
controller interpolation period, the simulation results 
demonstrate good performance of interpolating LPV 
control as shown in Figure 4. 

Concluding Remarks 

In this paper, an LPV control interpolation algo- 
rithm was proposed to achieve high-performance of 
gain-scheduled control. The proposed LPV control 
approach unified the systematic LPV control theory 
with interpolation technique. The stability of the 
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Time (s) 


(a) rotor displacement Xi , x 2 



Time (s) 



(b) control force Ui,u -2 


Pig. 3 Fixed parameter unit-step response with 
LPV controller. 

interpolated LPV control was achieved by enforcing 
continuity of Lyapunov function over entire parameter 
space. With modest increase of computational effort, 
the proposed LPV interpolation scheme improved the 
LPV control performance considerably. The newly de- 
veloped LPV control synthesis method was applied to 
the magnetic bearing control problem to reject un- 
balancing sinusoidal disturbances and accommodate 
changing rotor speed. Promising simulation results 
were obtained. 

For future research, an F-16 nonlinear model has 
been acquired from the NASA Langley Research Cen- 
ter and the proposed LPV control technique is cur- 
rently being applied to a highly maneuverable aircraft 
to improve its performance over the expanded flight 
envelope into nonlinear unsteady flight regimes. Be- 
cause aerodynamic controls may be insufficient for 
control at post-stall conditions, we would like to aug- 
ment the conventional aerodynamic surface controls 
using thrust vectoring (TV) control effectors. The 
LPV interpolated control design approach is capable 
of unifying the aerodynamic force and thrust force 
control law development in a unified framework, and 



(a) Rotor speed profile 



(b) Displacement: ap -solid, x 2 ”dash 


Fig. 4 Performance of interpolating LPV control 
for time-varying rotor speed. 

of providing guaranteed stability and performance for 
large flight envelope. 
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